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Large Hubbard U limit of the Kane-Mele model on a zigzag ribbon of honeycomb lattice near 
half-filling is studied via a renormalized mean-field theory. The ground state exhibits time-reversal 
symmetry (TRS) breaking d^ 2 _y 2 + idxy-waxe superconductivity. At large spin-orbit coupling, the 
Z 2 phase with non-trivial spin Chern number in the pure Kane-Mele model is persistent into the TRS 
broken state (called “spin-Chern phase”), and has two pairs of counter-propagating helical Majorana 
modes at the edges. As the spin-orbit coupling is reduced, the system undergoes a topological 
quantum phase transition from the spin-Chern to chiral superconducting states. Possible relevance 
of our results to adatom-doped graphene and irridate compounds is discussed. 

PACS numbers: 72.15.Qm, 7.23.-b, 03.65.Yz 


Introduction. Searching for topological states of quan¬ 
tum matters constitutes one of the central and funda¬ 
mental issues in condensed matter systems. The grow¬ 
ing interest in topological insulators (TIs), which sup¬ 
port gapless edge (or surface) states protected by time- 
reversal symmetry (TRS) while the bulk remains insulat¬ 
ing [1, 2], is one prime example. Of particular interest are 
topological superconductors which support gapless self¬ 
conjugate, charge-neutral fermionic quasi-particle excita¬ 
tions [3]. These excitations which reflect non-trivial topo¬ 
logical bulk properties are localized at the edges, known 
as Majorana fermions (MFs). 

Much effort has been put in searching for signatures 
of Majorana fermions in solid state materials. One¬ 
dimensional semiconductor nano-wires with strong SO 
coupling under a magnetic field proximity to a s-wave 
superconductor have been proposed theoretically to host 
MF at both ends of the wire [4-6], and also studied ex¬ 
perimentally [7-11]. Similar ideas have been proposed in 
2D systems where chiral MFs exist at the edges of spin- 
triplet, p—wave (odd-parity) superconductors [12-14]. 

While realization of the above systems relies on TRS 
breaking by the Zeeman field, time-reversal invariant 
topological superconductors (TRITOPs) [15-18] has re¬ 
cently been proposed to host two time-reversal pairs of 
helical MFs at edges in repulsively interacting SO cou¬ 
pled nano-wire proximity to either a s-wave [6, 16] or 
a d—wave [19] superconductor at each end of the wire. 
Proposals to realize TRITOPs in 2D systems include the 
spin-triplet Px i Wy superconductors [15], the bi-layer 
Rashba system [20], and in exciton condensates [21]. 

In this letter, we suggest a novel mechanism for re¬ 
alizing helical Majorana fermions in 2D topological su¬ 
perconductors with spin-singlet and TRS-breaking pair¬ 
ing gap- by directly doping correlated 2D quantum spin 


Hall insulators (QSHIs or 2D TIs) [22, 23]. A paradig¬ 
matic model for QSHIs is the Kane-Mele (KM) model 
[23, 24], which shows a non-trivial Z 2 topological (or 
spin Chern) number and supports helical topological 
edge states protected by TRS. The half-filled KM model 
with strong electron correlations is in the Mott-insulating 
(MI) phase [25], while superconductivity appears upon 
doping. Attractive candidates to realize correlated KM 
model include: graphene with enhanced SO coupling 
(^ 20meK) by doping with heavy adatoms, such as in¬ 
dium or thallium [26], Iridium-based honeycomb com¬ 
pounds X 2 lrOs{X = Na or Li) with strong SO coupling 
and electron correlations [25, 27]. 


Via renormalized mean-field theory (RMFT) approach 

[28] , we find the spin-singlet TRS breaking dx^-y ‘2 + idxy- 
wave superconductivity to appear at the ground state 
where the chiral edge states have been shown to occur 

[29] . Surprisingly, for sufficiently large SO coupling, in¬ 
stead of chiral edge states, we find zero-energy helical 
MFs to appear at each ribbon edge. We further show 
that the system exhibits a non-trivial Z 2 topological in¬ 
variant, which supports helical MFs at edges despite the 
presence of TRS breaking superconductivity. This seem¬ 
ingly un-expected feature comes as a result of persistence 
of spin-Chern phase in the pure Kane-Mele model in the 
superconducting state at finite dopings. The competi¬ 
tion between SO coupling and chiral d—wave supercon¬ 
ductivity leads to a novel spin-Chern to chiral topological 
quantum phase transition. 


Model Hamiltonian. The Hamiltonian of the Kane- 
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FIG. 1: (Color online) (a). Honeycomb lattice of a finite¬ 
sized zigzag ribbon of the tight-binding Kane-Mele t-J model 
with the ribbon size N = S being twice the number of zigzag 
chains along x—axis. Nearest-neighbor and next-nearest- 
neighbor lattice vectors are ^a=i, 2 , 3 , ai=i ,2 with an unit length 
of a, a, respectively. We set a = 1 here. The gray shaded re¬ 
gion represents for the super-unit-cell of the zigzag ribbon, 
which repeats itself along a:—axis. The filled (open) dots 
stand for the sites on sub-lattice A{B). The three phases for 
d+ id' pairing gap are dehned (see shaded green triangle) as: 
00 = 0,01(2) = — (+)27r/3. (b). Mean-field free energy Fmf 
versus doping S of doped Kane-Mele model on 2D periodic 
lattice for extended s-wave and d + zd'-wave superconducting 
pairing symmetries, (c). Mean-field variables Ao (strength of 
pairing gap) and X"field versus doping for (b). The parame¬ 
ters used in (b) and (c) are: J/t = 0.1, tsolt = 0.8. 


Mele t-J (KM-tJ) model is given by [23]: 

H = Hkm + Hj, 

Hkm = -tgt 

{ij),a i,a 

. A SO z a . 1 

+ 2 ^ l^ijC\ (7 Cj Fh.C. 

Hj = Jg,Y.^Si-Si-\nin^) (.1) 

(hi) 

where a stands for the spin index, (00) and 

refer to the nearest-neighbor (NN) and next-nearest- 
neighbor (NNN) sites, respectively (see Fig. 1 (a)). Here, 
Uij = 1 for 00 G A and Uij = —1 for i,0 G B in 
the SO coupling; Si refers to the electron spin opera¬ 
tor on site 0 defined as: Si = 1/2 ^ , 

rii = is the electron density operator, and the 

anti-ferromagnetic spin-exchange couping J ^ can 


be derived via the second-order perturbation from the 
Kane-Mele Hubbard model in the limit of a strong on¬ 
site Coulomb repulsion U ^ 030]. 

The Hj term has been known to favor the spin- 
singlet pairing. To address superconductivity of the 
model, we apply RMFT based on Gutzwiller projected 
single-occupancy constraint due to the proximity of the 
Mott insulating ground states, known to well describe 
the ground state of d—wave cuprate superconductors 
[28]. As a result, the hopping t term acquires a re¬ 
duction factor gt at a finite doping 5\ t ^ tgt with 
gt = 2J/(1 + J), while the spin-exchange J term gets 
enhances by a factor gs'. J gsJ with gs = 4/(1 + J)^. 
The spin-exchange J term within RMFT reads: Hj = 


^ + Hconst^ ^0, 




Ha = 




* + <5a 


* i+&J 


* + <5a 


h.c., 


+ h.c., He. 


Iv- F lA- F 

Hs + FGt] “ where a = 1,2,3, N, is 

the total number of sites, Xs ~ 

Based on the C^y symmetry of the underlying lattice, 
the pairing symmetry of 2 3 take the following 

forms [31]: (i) extended 5 —wave: A^ = At = Aq, (ii) 
dx‘ 2 -y ‘2 + idxy—vj 8 iYe (denoted also as d + id'): A^ = 
^ Aoe'^“-i with 0o = 0, 0i(2) = -(+)27r/3 (see 
Fig. 1) [32]. The Fourier transformed pairing gap Aj^ for 
a periodic 2D lattice reads: A^ = X]a=i 2 3 

The mean-field Hamiltonian 
on a periodic lattice in the basis of 0/. = 

(J J A A At ^tt At \T • 

V^A,/c’ '-A-k^ '-B-k^ '-A-k^ ^B,-k) 

given by the 8x8 matrix Aik’ 


Mk 

ht 

Ik 

ek 
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±7fc — fi Ck \ X _ f ^ Ak\ 

4 T7fc-hy’ A-fc 

/so [-sin(fcy) + 2cos{VSkx/‘2) sin(fcy/2)], 

itgt + x) 


a=l,2,3 


(• 2 ) 


The Hamiltonian Eq. (.2) possesses both the Particle- 
hole (PH) symmetry: C~^M.kC = —Ai-k^ C = r^K 
(with being the Pauli matrix on particle-hole space 
where K stands for complex conjugation as well as sub¬ 
lattice symmetry : Adk A4-k for CA,k (^B,k [31]- 
The matrix describing the KM model, shows TRS: 
T~^hkT = h-k where T = icr'^K is the time-reversal 
operator taking (c^,c^) to However, Adk 

breaks the TRS for dAid! superconducting order param- 
eter: = cos(7r/3)Ad^2_^2 (fc) + / sin(7r/3)Ad^^ (fc) 

with Arf^ 2 _^ 2 (fc) = 2Ao(e“®^«'/'^ — e*^i'/^cos(/i:3;/2-\/3)), 
^dxyik) = -2iAoe^'^y/‘^sm{kx/2\/3), and 
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T T 7^ ^d+^d mean-field free en¬ 
ergy reads: Fmf = “1^ S/c lii[cosh(^)] + with 

> 0 being positive eigenvalues and Ng the number of 
sites. We diagonalized the mean-field Hamiltonian Hk 
on a finite-sized zigzag ribbon with Ns = N/2 zigzag 
chains and N = 56 is set as the total number of sites 
along axis throughout the paper. 

Results. The mean-field variables are solved self- 
consistently by minimizing the free energy both for a 
periodic lattice and a finite-sized ribbon. On a 2D peri¬ 
odic lattice, the results as a function of doping are shown 
in Fig. 1(b) and (c). Compared to the TRS extended 
5—wave, we find dx 2 _y 2 ^idxy—wdive pairing is the ground 
state [32]. Same pairing symmetry has been reported in 
superconducting phase of the doped graphene in the ab¬ 
sence of the spin-orbit coupling [29, 31, 33, 34], which 
was argued to support two co-propagating chiral edge 
states at low energies with a non-trivial topological wind¬ 
ing number Ntknn = [35]. The superconducting 

transition temperature Tc is estimated as Tc ^ gt^o- 

On a finite-sized zigzag ribbon and at a generic dop¬ 
ing, the Bogoliubov quasi-particle dispersion shows four 
doubly-degenerate bulk bands (due to the Sz symme¬ 
try of our Hamiltonian) grouped in two pairs (see Fig. 
2(a)); it satisfies the particle-hole symmetry with 27r 
periodicity: E{kx — tt) = —E{kx + tt). At low dop¬ 
ings, the normal state Fermi surfaces enclose the Dirac 
points K± = (;^,±^) (see Fig. 2 (b)); the d + id' 
pairing strength is weak near K_ [32]. Surprisingly, 
in the regime of a strong SO coupling and weak pair- 
ing < \/3iso), we find the low energy excita¬ 

tions of our model support helical MFs at edges instead 
of chiral edge states as expected for a chiral d—wave 
superconductor. On a finite-sized zigzag ribbon, we 
find two Dirac-dispersed lines intersecting at momenta 
^ jZ where the Bogoliubov quasi-particle 

excitation energy vanishes, E{k^^) = 0 (see Fig. 2(a) 
and Fig. 3(a)) [36]. Near each of these gapless points, 
two pairs of two-fold degenerate states are generated via 
intersecting the two Dirac lines by a constant energy at 
two momentum points denoted as: 

These four degenerate states are located at the same 
edge. However, and are counter-propagating, 

while and are co-propagating (see Fig. 

3(b)). In the basis of T = [H^cJ, H^cj^, H^cJ^, H^ct'*^] 
with i = 1 • • • A/", the square magnitudes of two pairs 
of degenerate eigenstate wave-functions associated with 
the same k^^^ defined as |T(i)p = = 

(|i4(i)p, |'d(i)P, |'5(^)P, exhibit an exponential 

decay from both edges into the bulk, where R/L refers 
to the right/left moving state, and u{i)^u{i)^v{i)^v{i) are 
the corresponding matrix elements. 

These features are clearly different from the co- 
propagating chiral edge states realized either in the 
chiral d-wave superconductivity in doped graphene or 



FIG. 2: (Color online) (a) The Bogoliubov dispersion E{kx) 
(in unit of t) of doped KM-t J model on a zigzag ribbon with 
N — 56 for J/t — 0.1, tsolt — 0.8, and 5 — 0.05. (b). The 
spin-up Fermi surfaces in the normal state of the Kane-Mele 
model on 2D periodic lattice for iso jt = 1 at various dopings. 
The spin-down Fermi surfaces are obtained via A+ ^ K-. 
(c) The 3D density plot for the effective intra-band supercon¬ 
ducting gap function |A^“| of the Kane-Mele model on 2D 
periodic lattice (see text and Ref. 37). 


by proximity of a 5 —wave superconductor to a quan¬ 
tum anamolous Hall insulators [38]. Instead, the edge 
states we find fit well to the helical MFs described by 
the linearly-dispersed Hamiltonian defined by the Bo¬ 
goliubov quasi-particle operators as: = 


/c,2 k,i k,i —k,i k,i —k,i^ 


it = 


vl-d ui - + uIjcF^ where k = k - k^^, 


k,i —k,i 


k,i k,i 


k = kx, 7^^ with a = L^R refers to the Bogoliubov 
quasi-particle destruction operator defined by the co¬ 
herence factors pppp corresponding to the 
right-moving quasi-particle with “pseudo-spin” r =t (i 
) , and the summation over repeated site index i = 
1,---A/' is implied; similarly for • The pair of 

the degenerate wavefunctions 

can be expressed as formed by the coher¬ 
ence factors: T^’^(i) = ), T^’f(i) = 


(—14^ ^) ; similarly for the other doublet 

(i). It is clear from Fig. 3 (b) that the edge 
states (as well as the Bogoliubov 

operators form pairs (see pink (blue) 

curve in Fig. 3(b) for |T^’^(i)p(|T^’f(i)p)). Further¬ 
more, these Bogoliubov operators with linear dispersion 
satisfy ^_j^{—E) = 7|(^) with k = k — tt via PH sym¬ 
metry (see top and bottom panels of Fig. 3(b)). Hence, 
they can be regarded as examples of helical MFs at edges 
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FIG. 3: (Color online) (a) Bogoliubov excitation spec¬ 

trum of Fig. 2(a) near zero energy, (b) The square magni¬ 
tudes of two pairs of degenerate eigen-vectors |T^/^’'^(i)|^ = 
{\u(i)\‘^,\u{i)\‘^,\v{i)\‘^,\v(i)\‘^) (see text) for a fixed eigenen- 
ergy E = — ±0.014t with i running (left to right) from 

i = 1 (top edge) to i = 56 (bottom edge), corresponding to 
the helical Majorana fermions. Physical parameters are the 
same as in (a), (c) Schematic plot of the helical edge states 
in (b) for E = E^ where same color in (b) and (c) refers to 
the same state. Here, ZGNR refers to the Kane-Mele zigzag 
nano-ribbon. 


[15]; the MF zero-modes occur at /c = (or k = 0) 

where helical MFs are protected 

by an additional P-H symmetry: 7-fe(-E) = 'rliE) (see 
Fig. 3(b)). Our seemingly unexpected results have roots 
in the competition between TRS SO couping and TRS 
breaking chiral d—wave superconductivity. We find the 
TRS protected Z 2 QSH insulating phase of the pure un¬ 
doped Kane-Mele model persists up to a finite doping and 
a finite pairing gap in the chiral d—wave superconducting 
phase. 

To gain more understanding, we compute the Z 2 topo¬ 
logical invariant (spin Chern number) Nw = ^{Cn — Cfi) 
with Cn = ^ where the integral is done 

in the first Brillouin zone (BZ), the field strength Fi 2 {k) 
and the associated Berry’s connection A{k) are defined 
as: Fi 2 {k) = diA 2 {k)-d 2 Ai{k) and = {n{k)\d^\n{k)) 
with \n{k)) being the normalized eigenvector of the nth 
band [39]. is defined similarly with n being the 

corresponding degenerate band associated with the n- 
th band via the transformation Sz In the strong 

SO coupling regime, tso ^ ^0 and for a sizable range 
of doping around 1/2-filling, the spin-Chern phase pre¬ 
vails and we find Cn = —Cn or Nw = il [37], same 
result as that for the Z 2 TRITOPs supporting helical 
MFs at edges [1, 15, 40, 41]. In this limit, our system is 
well approximated by the effective spin-singlet Px ± ipy 


(a) Chiral spin-Chern Chiral 
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FIG. 4: (Golor online) (a) Topological phase diagram of 
doped Kane-Mele t-J model on 2D periodic honeycomb lattice 
as a function of p for fixed Ao/t = 0.3, tso It = 0.9. Here, 
we set X = 0, = ps = 1 for simplicity, and Epc ^ ±V^tso 

refers to the critical chemical potential at the spin-Ghern- 
chiral phase transition. The values of Ao and p are tuned in 
a non-self-consistent way. (b) Energy dispersion of the two 
bulk bands close to zero energy. The bulk band gap closes 
at the phase boundary p ~ —pc at one Dirac point (case (ii) 
with pjt — —1.409, S = 0.43), while it remains hnite on either 
of the two phases (cases (i) with p/t = —1.8, (5 = 0.58 and 
case (hi) with p/t = —1.2, (5 = 0.34). 


superconductivity near the two Dirac points K±. This 
can be seen when re-expressing the superconducting pair¬ 
ings in terms of the electron operators '0±,/c which diag- 
onaliz the tight-binding KM Hamiltonian [32, 37]. In 
this basis, the intra-band pairing dom¬ 

inates at ground state (see Fig. 2 (c) and Ref. 37). 
Near K± points with q± = K± + (±^x + we find 
A“^“ ^ ±Ao(qy ± iqx), resembling the case of a TRI¬ 
TOPs. In the opposite limit for tso ^ Aq or suffi¬ 
ciently large doping where the chiral d-wave pairing dom¬ 
inates, however, we recover the chiral superconductivity: 
Nw = 0 and Cn = Cn = 1, equivalent to the case of 
doped graphene [29, 31]. A novel spin-Chern-to-chiral 
topological quantum phase transition is found as A^/tso 
or p/tso is varied [42]. The generic phase diagram by 
tuning p (in a non-self-consistent way) at a fixed Ao/tso 
is shown in Fig.4 (a). For Aq ^ tso^ we find the critical 
values of p being at p = Epc r\j E^/3tsO‘ The bulk band 
gap closes only at the phase transition, while it remains 
open in either phase [42-44] (see Fig. 4(b)). The similar 
persistence of spin-Chern phase in a TRS breaking mag¬ 
netic field has been observed experimentally in Ref. 45. 
The spin-Chern phase of our system belongs to class D 
topological superconductors, distinct from the TRITOPs 
[19, 46, 47]. 

Conclusions. We demonstrate here for the first 
time a 2D spin-singlet topological superconductor with 
non-trivial spin Chern number in doped correlated Kane- 
Mele model, which supports helical counter-propagating 
Majorana zero modes despite the d^id' superconducting 
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pairing gap breaks TRS. This seemingly unexpected fea¬ 
ture comes as a result of persistence of spin-Chern phase 
of the pure Kane-Mele model in the superconducting 
state upon doping. As T ^ 0, distinct differential con¬ 
ductance spectrum for each pair of Majorana zero mode 
through differential Andreev conductance in the normal- 
metal/superconducting (NS) junction is expected[29, 31]. 
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In the Supplementary Materials, we provide some of the details in the main text. 


PACS numbers: 72.15.Qm, 7.23.-b, 03.65.Yz 


I. THE HAMILTONIAN OF THE KANE-MELE ing matrix form: 
T-J MODEL ON A ZIGZAG RIBBON 


Here, we provide in details the Hamiltonian of the 
Kane-Mele t-J model on zigzag ribbon. The Hamiltonian 
matrix is AN x AN in size with A/2 being the number 
of zigzag chains along x-axis. Each zigzag chain passes 
through either A sites or B sites. The ket state of the 
Hamiltonian is 


1^) — * * * Cat^) 

(I.l) 


for i = 3,5,7---A-3 

2 k 

Hso^cFcrij 2 , i) — Hi 2 ), 

2 k 

Hso,aa{U i “h 2) — (J — tgoSifli^ — ^, 

Hso,(j(j{j Y 2, i) = HSQ (jcjif-) ^ 2), 

T 

2 


The Hamiltonian comprises three terms 


H = HoPHgoPHA. (L2) 


Here, Hq represents the nearest-neighbor hopping terms 
in the matrix form of (* * * • • • )H{’ • • Cja- • • • )^, given 

by: 


for i = 4, 6, 8 • • • A - 2 

2 k 

Hgo^o-crijC 2) (J — — ) , 

Hgo^crcrij 2, i) = HgQ ucriii i 2), 

2 k 

Hgo,aa{U i Y 2) — (J — tgoSifll^ — ^ , 

Hgo^<j(j{j Y 2, i) = HgQ_^(j(j{i, i 2), 

2 . 

Hgo^<j(j{j 1 ^'^tgoSifii^k^, 

o 

and the other matrix elements read 


for i = 1,3, 5 • • • A — 1, 

k 

-^0,crcr(^5 i Y 1) 2tc05( —), 

-^0,crcr(^ Y lA) — i “h 1) 

for i = 2,4,6---A, 

HQ^(j(j{i,i Y 1) — t, 

-^0,crcr(^ Y lA) — i “b 1), 


(1.3) 


-^SOjCrcr (I 5 1 ) 2 5 

2 

-^SOjCrcr (2, 2) (J — tgQSifli^k^ , 

2 

Hso,aa 1 ) (7 -H gQsiniJCj, 

2 

Hgo,aa{^ 1 


(1.4) 


where a =t,i (Y or —), and k = kx- The matrix el¬ 
ements of the form (• • • • • • )^(’ ’' related 

to (• • • • • • )H{' • • Cja- • • • )^ by a minus sign. The Hgo 

describes the intrinsic spin-orbit couping with the follow- 


where the matrix elements with the abnormal order is mi¬ 
nus the normal order elements. The Ha is the supercon¬ 
ducting pairing gap with the matrix elements of the form 

or 
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FIG. 1: (Color online) (a) Superconducting gap as a function of position in the zigzag ribbon for iV = 26 (or N/2 = 13 

zigzag chains), (b) Square magnitude of wave functions |^(i)|^ at the edge in the superconducting states with an energy 
Ejt — —0.613 as a function of position i — 1 • • • N in the zigzag ribbon with N = 56. Here, we take the uniform mean-field 
solution for (i) = 0.22 for N = 56. The parameters used in (a) and (b) are: J/t = 0.1, 6 = 0.15, tso/t = 0.5. 


given by: 

for i = 1,3, 5 • • • — 1 

HA{i,i-lE3N) = Ao 
— 1g + 3N) = Aq 
Ha{'^^ i + 1 + 3N) = Ai + A 2 
Ha{'^ + 1, i + 3N) = Ai + A 2 
similarly 

Ha{'^ T a/”, i — 1 2N) = — Aq 

Ha^'i ~ 1 + i + 2N) = — Aq 

HA{i + A, i + 1 + 2A) = -Ai - A 2 

Ha{'^ + 1 + A, i + 2A) = —Ai — A 2 

HA{i + 2A, i - 1 + A) = HA{i + A, i - 1 + 2A)* 

HA{i - 1 + 2A, i + A) = HA{i - 1 + A, i + 2A)* 

HA{i + 2A, i + 1 + A) = HA{i + A, i + 1 + 2A)* 

HA{i + 1 + 2A, i + A) = HA{i + 1 + A, i + 2A)* 

HA{i + 3A,i - 1) = AA(iG - 1 + 3A)* 

HA{i - 1 + 3A,i) = HA{i - 1, i + 3A)* 

Aa(^ + 3A, i + 1) = Aa(^, i + 1 + 3A)* 

AA(i + 1 + 3A, i) = AA(i + l,i + 3A)*, 
and the other matrix elements read 
Aa(1,2 + 3A) = A 1 +A 2 
Aa(2, 1 + 3A) = Aa(1, 2 + 3A) 

Aa(1 + a, 2 + 2A) = -Aa(1, 2 + 3A) 

Aa(2 + a, 1 + 2A) = -Aa(2, 1 + 3A) 

Aa(1 + 2A, 2 + A) = Aa(1 + a, 2 + 2A)* 

Aa(2 + 2A, 1 + A) = Aa(2 + a, 1 + 2A)* 

Aa(1 + 3A,2) = Aa(1,2 + 3A)* 

Aa(2 + 3A, 1) = Aa(2, 1 + 3A)* 


where Aq is the pairing order {ci^aCj^-a)^ Aq, Ai and A 2 
correspond to the pairing order parameters carrying dif¬ 
ferent phases 0,^ and ^ in d+zd'-wave superconducting 
states, respectively. 


II. ADDITIONAL RESULTS OF THE 
KANE-MELE T-J MODEL ON A ZIGZAG 
RIBBON 

Here, we provide additional results based on our nu¬ 
merical calculations on the doped Kane-Mele t — J model 
via RMFT. 


A. Superconducting states at the edges 

In the Bogoliubov quasi-particle spectrum on the 
zigzag ribbon (see Fig. 2 in the main text), we find 
two linear-dispersed spectrum crossed at kx ^ tt with 
the energy E = EEq {Eq > 0). They are closely related 
to the edge states of the pure Kane-Mele model. The 
negative energy state {E < 0) corresponds to the super¬ 
conducting state at the edges as it supports finite values 
in superconducting order parameter Aq at the edges and 
their corresponding eigen-vectors are distributed mainly 
at the edges (see Fig. 1 (b)). 

To investigate further the position dependence of the 
superconducting gap , we generalize our mean- 

field calculations to allow for spatially varying supercon¬ 
ducting gap A{i) (instead of the gap with an uniform 
magnitude). As shown in Fig. 1 (a), the magnitude of 
gap is enhanced at edges and stay at a constant value 
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FIG. 2: Bogoliubov excitation spectrum of the Kane-Mele t-J model for different dopings. The parameters used here are: 
J/t — 0.1, tsojt = 0.5. 


in the bulk. This can be understood as the electronic 
density of states (DOS) of the un-doped pure KM ribbon 
is enhanced at the edges due to the presence of the heli¬ 
cal edge states. As a result, the superconducting pairing 
strength is enhanced at edges. However, the inhomogene- 
ity of superconducting gap at the edges concerns only 
with the bulk states at a sizable negative energy. Since 
we focus in this work on the Majorana fermions near 
zero energy, far away from the superconducting states at 
edges, we therefore ignore this spatial inhomogeneity and 
consider only the superconducting gap with an uniform 
magnitude within our RMFT approach. 


B. Doping evolution of the Bogoliubov excitation 
spectrum 

Here, we provide the doping evolution of the Bogoli¬ 
ubov excitation spectrum as shown in Fig. 2. The 
spin-Chern phase with helical Majorana modes survive 
at low dopings. For J/t = 0.1,t5oA = survives 

for 0 < 6 < 0.2. Above 6 = 0.2, the superconducting 
gap vanishes and the system is in the normal state. On 
the other hand, the chiral phase will prevail in the op¬ 
posite regime tso ^ for fow dopings. However, the 
chiral superconductivity has been addressed extensively 
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(see Refs. 29, 35 in the main text). We will leave this 
issue out of the focus of our present work. 


III. Z2 TOPOLOGICAL INVARIANT 


In this section, we provide details on obtaining the Z 2 
topological number (or the spin Chern number) in the 
doped Kane-Mele t-J model following the approach in 
Ref. 2. First, the TKNN number, a topological (Chern) 
number, for the n—th band is defined as^: 

Cn = 7^ [ de Vk X A(k), (III.l) 

where A(k) = (n(k)|Vk|^(k)), |n(k)) is the Bloch state 
in the n-th band up to a normalization coefficient, 
which satisfies the Schroedinger equation i^(k)|n(k)) = 
F^^(k)|n(k)). The existence of A(k) implies that there is 
a non-trivial magnetic field across the first Brillouin Zone 
(FBZ). The FBZ of the honeycomb lattice we choose for 
the integration Eq. (III.l) is illustrated in Fig. 3, the 
parallelogram (area enclosed by red lines) spanned by 
vectors qi and q 2 in momentum space. 

In order to numerically calculate the TKNN number 
covered by the entire FBZ (area enclosed by the solid 
red lines) in Fig. 3, we first discretize it into N — 1 x 
N — 1 small flakes of parallelograms such as the smaller 
parallelogram abed (enclosed by solid black lines) in Fig. 
3. The position of a discretized lattice point (that is the 
vertices of the small parallelograms ) can be specified by 
a vector k^j = (, kj 2 ), for (i^j) G [1 , V] with its 
components 

kil = ~ [("^^ ~ 1) ^ Qlx {j ~ f) ^ Q2x\ 5 (III.2) 

^ ^ «2y] , (111.3) 


where Qix and Qiy are the components of the vectors qi 
and q 2 with unit lattice spacing |a| = 1 , they are given 
by 


qi = {qix 



q2 = {q2x , q2y) 


27r 2-\/37r\ 

T ’ ) ■ 


(IIL4) 

(III.5) 


The phase difference A0i(ky) of (k) = (k|n(k)) 
travelling from a lattice point k^^j to can be com¬ 

puted to be 


A6>i(kij) = — iln 


(kij) (ki+ij) I 


= -i\nUi (ky), 


where 


Ui (k,,,) 


(km,j) 


(IIL6) 


(111.7) 


The subscript ”1” of A^i and Ui (k^j) represents that 
4^^ travels along the direction of qi. Sum over AO of the 
four sides of the small parallelogram abed in Fig. 3 can 
one obtain the total phase shift of the Bloch wavefunc- 
tion on the n-th band after circling around that small 
parallelogram. Compare with Eq. (III.l), we define the 
lattice 1/(1) gauge field strength F-j by 

F" = Inf/i (ki.,) U2 (k,+i,,) I/f 1 (k,,,+i) f/2-i (k,,,) 

(III.8) 

Einally, we can calculate the TKNN number on the lat¬ 
tice associated the n-th band by 



a = (*i/ . *2y) 

* = (*i,>i . *2y) 
c = (^,+i , ^2y+i) 
d = {ku , 


FIG. 3: The parallelogram enclosed by red solid lines is the 
area of entire FBZ with F = (kx — 0, ky = 0). The small flake 
of parallelogram abed is the smallest unit after discretizing 
the FBZ for numerical calculation. 


= (nw) 

For the same band we can further define a topological 
number, so-called the spin Chern number, as: 

Nw = (III, 10) 

where Cn is the TKNN number of the other degenerate 
eigenvector of the same band. Numerically, we obtain 
Cfi = — Cn, thus Nw is either 1 or —1. We checked 
numerically that for the parameters used in Fig. 3 of the 
main text \Nw\ approaches to 1. 

In the other extreme limit, A tso^ we found that 
Ntknn = Cn F Cfi = ±2, and Nw = 0, as expected for 
d + id' superconductivity in doped graphene^. Further 
studies suggest the existence of a quantum phase transi¬ 
tion at a critical ratio of A/tso separating the spin-Chern 
phase with Nw = ±1 from the chiral (TKNN) phase with 
Ntknn — (see Ref. 3). 


IV. EFFECTIVE TRS SUPERCONDUCTIVITY 
NEAR THE DIRAC POINTS 

In this section, we derive the low-energy effective TRS 
superconducting pairing of our model near Dirac points. 
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We re-express the d + id' superconducting order param¬ 
eter in terms of the electron operators '0±,/c which di¬ 
agonalizes the tight-binding KM Hamiltonian with the 
corresponding eigenvalue -h 7 ^ — yu) (see 

Ref. 5): 


^A,k _ 

j ■ 


Pk Pk 


and 


4,/c 

^B.k 


Ak PP 


<,k 

pP-,k 


Ek 

V’i.fe 


(IV.l) 


(IV.2) 


where 

_Nd_ 


= P 


± ^h(lk±^Jel+^l) 

and {a^y = -cPp^, P^ = 


P 


± _ 


v/|«fcP + (7fe±A/«fc+7fc)^’ 

The d + id' superconducting pairing term can be 
re-written as: 

i^A = k ~ ^^Ak'^^Bk) + 

k 

= EAp~y-,k^Ek + ^k^y+,k^^Ek 


+ ^-AEp^Ek + ^p^AEkyi\k\ + h.c; 


. ++(—) 


^k l^/c ll^/c 


^+,/c 


^ ^^d+id'^d+id' 

A+”(”+^ = - At+y‘^'], 


(IV.3) 


where with 0 {a^) being th ephase 

of a^. In the presence of a finite SO coupling 
and hole-doping, the intra-band pairing of the lower 
band with pairing gap function domi¬ 
nates. The gap function near Dirac points K± = 

with q± = K± + (Eqx^qy) behaves to lead¬ 
ing order in |g| as an effective spin-singlet Px ± ip^-like 
superconductivity: 

Ag^ ^ —aQA.Qi[qx^iqy) (IV.4) 


with ao = ^ —. Hence, our system behaves as an 

12 v Stso 

effective TRITOPs. In deriving Eq. (IV.4), we have used 
the following results: 


- ^ ^q- 

^ 2^/3^5o’ 

aq_ « e^-Zk?-!- (IV.5) 

We have also used the approximated forms of d id' 
pairing gap (g)^ near the Dirac points, shown to 

be effectively a mixture of 5 —wave like (near Kj^) and 
Px + ipy—wave like (near K_) superconducting pairing^: 

A^+^^'(g+) « 3Aoe^^^/^ 

A‘'+*‘''(</_) « ^-Aoe^^/\qy-iqy. (IV.6) 
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